Modern compact star observations and the 
quark matter equation of state 



> 
o 

o 



T. Klahn, a ' b ' c D. Blaschke, b ' d ' e F. Sandin, f Ch. Fuchs,§ 
A. Faessler, g H. Grigorian, h;1 G. R6pke, b J. TriimperJ 

Physics Division, Argonne National Laboratory, Argonne, IL 60439-4843, USA 



o 

o, 

CN , b Institut fur Physik, Universitat Rostock, D- 18051 Rostock, Germany 



c Gesellschaft fur Schwerionenforschung mbH ( GSI), 
D- 64291 Darmstadt, Germany 

d Institute for Theoretical Physics, University of Wroclaw 



PL-50-204 Wroclaw, Poland 

e Bogoliubov Laboratory of Theoretical Physics, 
£> ' Joint Institute for Nuclear Research, RU-141980 Dubna, Russia 

t> ! f 

Department of Physics, Lulea University of Technology, SE-97187 Lulea, Sweden 

. g Institut fur Theoretische Physik, Universitat Tubingen, 

<^> D-72076 Tubingen, Germany 

\ h Laboratory for Information Technologies, 

Joint Institute for Nuclear Research, RU-141980 Dubna, Russia 

1 Department of Physics, Yerevan State University, 375049 Yerevan, Armenia 

O I J Ma.-Planck-InsUtut fur E X traterrestnsche Phys % k, D-857 4 l GarcUng, Germany 

> : 

^ , Abstract 

H . 

We present a hybrid equation of state (EoS) for dense matter that satisfies phe- 
nomenological constraints from modern compact star (CS) observations which in- 
dicate high maximum masses (M ~ 2M ) and large radii (R > 12 km). The corre- 
sponding isospin symmetric EoS is consistent with flow data analyses of heavy-ion 
collisions and a deconfinement transition at ~ 0.55 fm~ 3 . The quark matter phase 
is described by a 3-flavor Nambu-Jona-Lasinio model that accounts for scalar di- 
quark condensation and vector meson interactions while the nuclear matter phase 
is obtained within the Dirac-Brueckner-Hartree-Fock (DBHF) approach using the 
Bonn-A potential. We demonstrate that both pure neutron stars and neutron stars 
with quark matter cores are consistent with modern CS observations. Hybrid star 
configurations with a CFL quark core are unstable within the present model. 
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1 Introduction 



Recently, new observational limits for the mass and the mass-radius relation- 
ship of compact stars have been obtained which provide stringent constraints 
on the equation of state of strongly interacting matter at high densities, see [T] 
and references therein. In particular, the high mass of M — 2.1 ± 0.2 M for 
the pulsar J0751+1807 in a neutron star - white dwarf binary system [2J and 
the large radius of R > 12 km for the isolated neutron star RX J1856. 5-3754 
(shorthand: RX J1856) [3] point to a stiff equation of state at high densities. 
Measurements of high masses are also reported for compact stars in low-mass 
X-ray binaries (LMXBs) as, e.g., M = 2.0 ±0.1 M Q for the compact object in 
4U 1636-536 @|. For another LMXB, EXO 0748-676, constraints for the mass 
M > 2.10 ± 0.28 M & and the radius R > 13.8 ± 0.18 km have been reported 
[5]. The status of these data is, however, unclear since the observation of a 
gravitational redshift z = 0.35 in the X-ray burst spectra |6j could not be 
confirmed thereafter despite numerous attempts [7]. We exclude possible con- 
straints from LMXBs from the discussion in the present paper as their status 
is not settled and they would not tighten the mass and mass-radius limits 
provided by J0751+1807 and RX J1856, respectively. It has been argued [3"f5] 
that deconfined quark matter cannot exist in the centers of compact stars with 
masses and radii as reported for these objects. In view of recent works on the 
quark matter EoS, however, this claim appears to be premature [8]. 

In the present paper, we demonstrate that the present-day knowledge of hy- 
drodynamical properties of dense matter allows to construct hybrid EsoS with 
a critical density of the deconfinement phase transition low enough to allow 
for extended quark matter cores and stiff enough to comply with the new mass 
measurements of compact stars. It has been shown recently by Alford et al. [H] 
that hybrid stars can masquerade as neutron stars once the parameters of a 
generic phenomenological quark matter EoS have been chosen appropriately. 
While in [9] the APR EoS [10] for the hadronic phase has been used, we base 
our investigation on a nuclear EoS obtained from the ab initio DBHF ap- 
proach using the Bonn A potential [TT]fT2"] which results in star configurations 
with larger radii and masses. The DBHF EoS is soft at moderate densities 
(compressibility K=230 MeV) [1T][1~3"] but tends to become stiffer at high den- 
sities. At densities up to 2-3 times nuclear saturation density it is in agreement 
with constraints from heavy ion collisions based on collective flow [T4"1[T5] and 
kaon production [16]. However, at higher densities this EoS seems to be too 
repulsive concernig the flow constraint. As we will show in the present paper, 
the problems of this EoS with an early onset of the nuclear direct Urea (DU) 
process and a violation of the flow constraint for heavy-ion collisions at high 
densities can be solved by adopting a phase transition to quark matter. 

We have no first principle information from QCD about the quark matter EoS 
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in the nonperturbative domain close to the chiral/ deconfinement transition at 
zero temperature and finite density which would be required for an ab-initio 
study of the problem whether deconfined quark matter can exist in neutron 
stars or not. Therefore it is desireable to develop microscopic approaches to 
the quark matter EoS on the basis of effective models implementing, as far 
as possible, QCD symmetries into the model Lagrangian. Such an approach 
should in particular be capable of describing the phenomena of dynamical 
chiral symmetry breaking/restoration and color superconductivity in dense 
quark matter [T7], going thus beyond the traditional bag model approaches 
widely used for phenomenological studies of quark matter in compact stars, see 
[l"8][l"9] and references therein. A powerful modeling of QCD in the perturbative 
and nonperturbative domain is obtained using Dyson-Schwinger equations 
[2~0"f2T] which recently have been extended also to the finite density domain 
[22f23f24j . They appear as most promising candidate theories for a future 
systematic QCD modeling of dense QCD matter in compact stars. 

For the present study, however, we will employ the Lagrangian of a Nambu- 
Jona-Lasinio (NJL) model with chiral symmetry which is dynamically broken 
in the nonperturbative vacuum and restored at finite temperatures in accor- 
dance with lattice QCD simulations. Therefore, the application of the NJL 
model to finite densities where presently no reliable lattice QCD simulations 
exist, can be regarded as a step beyond the bag model phenomenology to- 
wards a dynamical approach to nonperturbative quark matter effects in com- 
pact stars, see [25] and references therein. Although the standard NJL model 
is quite simple, based on a few parameters, such as the momentum cutoff 
and coupling constants for a set of interaction channels, it offers possibili- 
ties for generalizations, like a density-dependence of the cutoff parameter [26] . 
momentum-dependent formfactors [27] and an infrared cutoff to account for 
confinement [28129] . 

In this work, we will use a standard three-flavor NJL model, where in contrast 
to the approach used in Ref. [9] we base our investigation on selfconsistently 
determined quark masses and pairing gaps [30], similar to the parallel devel- 
opments in Refs. [3"Tf 321133] . This approach has the advantage that it allows, 
e.g., to distinguish two- and three-flavor phases in quark matter (for a first dis- 
cussion, see [34J) and conclusions about the presence of gapless phases at zero 
temperature as a function of the coupling strengths in the current-current- 
type interaction of the model Lagrangian [35] . Moreover, we will investigate 
the question of the stability of neutron stars with a color superconducting 
quark matter core in the celebrated CFL phase, for which a number of phe- 
nomenological applications have been studied, in particular the cooling prob- 
lem [361137113 8 39 j. gamma-ray bursts [40|41] . and superbursts [42] . We will 
confirm in this work the earlier result by Baldo et al. [13] and Buballa et al. 
[44J that a CFL quark matter core renders the hybrid star unstable. 
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Here we generalize the model [30] by including an isoscalar vector meson 
current which, similar to the Walecka model for nuclear matter, leads to a 
stiffening of the quark matter EoS. Increasing the scalar diquark coupling 
constant leads to a lowering of the phase transition density. It is the aim 
of the present work to determine the unknown coupling strengths in both 
these channels such that an optimal hybrid star EoS is obtained. It fulfills all 
recently developed constraints from modern compact star observations [Tf5] 
while providing sufficient softness of the isospin-symmetric limit of this EoS 
as required from the analysis of heavy-ion collision transverse flow data [TUTo] 
and K + production data [IB] . 



2 Equation of state 



The thermodynamics of the deconfined quark matter phase is described within 
a three-flavor quark model of Nambu-Jona-Lasinio (NJL) type. The path- 
integral representation of the partition function is given by 



Z{T,fi) = f VqVqexp I ( dr j d 3 x [q (i$-m + /t7°) q + C 
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(2) 



where fx and rh = diagj(m M , m^, m s ) are the diagonal chemical potential and 
current quark mass matrices. For a = 0, r = (2/3) 1/,2 l/, otherwise r a and 
A a are Gell-Mann matrices acting in flavor and color spaces, respectively. 
C = 27 2 7° is the charge conjugation operator and q = q^°. G$, rjv, and 
t)d determine the coupling strengths of the interactions. 



The interaction terms represent current-current interactions in the color sin- 
glet scalar and vector meson channels, and the scalar color and flavor an- 
titriplet diquark channel. In the choice of the four-fermion interaction chan- 
nels we have omitted the chiral partner terms which would be necessary to 
establish the chiral symmetry of the Lagrangian as their contributions to the 
thermodynamical potential vanish at the mean-field (Hartree) level to which 
we restrict ourselves in the present paper. The model is similar to the models 
in [30ll31f32j . except that we include also the isoscalar vector meson channel. 
We follow the argument given in [30], that the Ua(1) symmetry breaking in 
the pseudoscalar isoscalar meson sector is dominated by quantum fluctuations 
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and no 't Hooft determinant interaction needs to be adopted for its realization. 

After bosonization using Hubbard- Stratonovich transformations, we obtain an 
exact transformation of the original partition function (JTJ). The transformed 
expression constitutes the starting point for powerful approximations, defined 
as truncations of the Taylor expanded action functional to different orders 
in the collective boson fields. In the following, we use the mean-field (MF) 
approximation. This means that the bosonic functional integrals are omitted 
and the collective fields are fixed at the extremum of the action. The corre- 
sponding mean-field thermodynamic potential, from which all thermodynamic 
quantities can be derived, is given by 



Q M f(T, /i) 



pv 



8Gc 



In Z MF (T, fj) 

i=u,d,s ^ V a=2,5,7 



" / (^5 t [ E ° + 2Tln (1 + + ft - n. . 



(3) 



Here, Qi is the thermodynamic potential for electrons and muons, and f2o is 
a divergent term that is subtracted in order to get zero pressure and energy 
density in vacuum (T = /i = 0). The quasiparticle dispersion relations, E a (p), 
are the eigenvalues of the Hermitean matrix 



M 



+ 7°/i* 7 5 r A A A A AA 
--f 5 T A X A A* AA -f T -p + m* - 7°/i* 



(4) 



in color, flavor, Dirac, and Nambu-Gorkov space. Here, are the diquark 
gaps, m* is the diagonal renormalized mass matrix and fi* the renormal- 
ized chemical potential matrix, fi* = diagy(/i u — Gsr/v^o, Hd — GsVv^o^s — 
GsVv4>o)- The gaps and the renormalized masses are determined by mini- 
mization of the mean-field thermodynamic potential ([3D- We have to obey 
constraints of charge neutrality which depend on the application we consider. 
In the (approximately) isospin symmetric situation of a heavy-ion collision, 
the color charges are neutralized, while the electric charge in general is non- 
zero. For matter in /^-equilibrium, also the electric charge is neutralized. For 
further details, sec [30.3L32J. 

We consider t]d as a free parameter of the quark matter model, to be tuned 
with the present phenomenological constraints on the high-density EoS. Sim- 
ilarly, the relation between the coupling in the scalar and vector meson chan- 
nels, rjv, is considered as a free parameter of the model. The remaining degrees 
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of freedom are fixed according to the NJL model parameterization in Table I 
of [13], where a fit to low-energy phenomenological results has been made. 

A consistent relativistic approach to the quark hadron phase transition where 
the hadrons appear as bound states of quarks is not yet developed. First steps 
in the direction of such a unified approach to quark-hadron matter have been 
accomplished within the NJL model in [29]. In this paper, however, the role 
of quark exchange interactions between hadrons (Pauli principle on the quark 
level) has yet been disregarded. As has been demonstrated within a nonrela- 
tivistic potential model approach, these contributions may be essential for a 
proper understanding of the short-range repulsion [16] as well as the asym- 
metry energy at high-densities [47J. In the present work we apply a so-called 
two-phase description where the nuclear matter phase is described within the 
relativistic Dirac-Brueckner-Hartree-Fock (DBHF) approach considered in [1] 
and the transition to the quark matter phase given above is obtained by a 
Maxwell construction. The critical chemical potential of the phase transition 
is obtained from the equality of hadronic and quark matter pressures. A dis- 
cussion of the reliability of the Maxwell construction for the case of a set of 
conserved charges is discussed in [H]. 

The baryon density as derivative of the pressure with respect to the bary- 
ochemical potential exhibits a jump at the phase transition, as shown for 
isospin-symmetric matter in Fig. [TJ As can be seen in that Figure, a slight 
variation of the quark matter model parameters t]d and rjv results in consid- 
erable changes of the critical density for the phase transition and the behaviour 
of the pressure (stiffness) at high densities. There appears the problem of a 
proper choice of these parameters which we suggest to solve by applying the 
flow constraint [H] to symmetric matter, shown as the hatched area in FigJTJ 
At first we fix the vector coupling by demanding that the high density behavior 
of the hybrid EoS should be as stiff as possible but still in accordance with the 
flow constraint. We obtain rjy = 0.50, independent of the choice of the scalar 
diquark coupling. The latter we want to determine such that the problem of 
the violation of the flow constraint for the DBHF EoS in symmetric nuclear 
matter at high densities is resolved by the phase transition to quark matter. 
The optimal choice for rjo is thus between 1.02 and 1.03. In the following we 
will investigate the compatibility of the now defined hybrid star equation of 
state with CS constraints. 



3 Astrophysical constraints on the high-density EoS 

Recently, observations of compact objects have resulted in new limits for 
masses and radii which put stringent constraints on the high-density behaviour 
of the nuclear matter EoS, see pQ. 



6 




0.2 0.4 0.6 0.8 1 

n [fm 3 ] 

Fig. 1. Pressure vs. density of the isospin symmetric EoS for different values of the 
relative coupling strengths r]£> and rjy. The behaviour of elliptic flow in heavy- ion 
collisions is related to the EoS of isospin symmetric matter. The upper and lower 
limits for the stiffness deduced from such analyses are indicated in the figure (shaded 
region). The quark matter EoS favored by the flow constraint has a vector coupling 
r]y = 0.50 and a diquark coupling between tjd = 1.02 (blue solid line) and rjo = 1-03 
(black fat solid line); results for four intermediate values rj£> = 1.022 . . . 1.028 are 
also shown (thin solid lines). 

Particularly demanding data come from the pulsar PSR J0751+1807 with a 
lower mass limit of M > 1.9 M [2], and the isolated compact star RX J1856 
with a mass-radius relationship supporting a radius exceeding 13.5 km for 
typical masses below 1.4 M or masses above 1.9 M for stars with radii 
R < 12 km PJ. In Fig[2] we display these constraints together with lines 
of constant gravitational redshift as, e.g., for the putative surface redshift 
measurement |6J of EXO 0748-676. 

The above constraints have to be explained by any reliable CS EoS, i.e. the 
mass radius relation resulting from a corresponding solution of the Tolman- 
Oppenheimer-Volkoff equations has to touch each of the regions shown in 
Fig. [2j This is well fulfilled for the purely hadronic DBHF EoS. 

It is widely assumed that if quark matter would exist in CSs, the maximum 
mass would be significantly lower than for nuclear matter stars (NMS). This 
argumentation has been used to claim that quark matter in neutron stars is 
in contradiction with observations [5]. 
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As we will show in this work, large hybrid star masses can be obtained for 
sufficiently stiff quark matter EsoS. In this case, the corresponding hybrid 
(NJL+DBHF) QCS sequence is not necessarily ruled out by phenomenology. 
The stiffness of the quark matter EoS can be significantly increased when the 
vector meson interaction of the NJL model introduced in Section [2] is active. 
The maximum value of the vector coupling which is still in accordance with 
the upper limit extrapolation of the flow constraint, r)v = 0.50, see Fig. 1, 
allows a maximum mass of 2.1 M Q . With this choice the constraints from 
PSR J0751+1807 and RX J1856 displayed in Fig. 2 can be fulfilled. 

The maximum mass is rather inert to changes of the diquark coupling t)d 
whereas the critical mass for the occurrence of a quark matter core gets sig- 
nificantly lowered by increasing the value of rj^,. For example, the choice of rju 
in the range 1.02 — 1.03 corresponds to critical star masses from 1.35 M & to 
1.0 M , see Figs. 2 and 4. 

We find for all hybrid EoS studied in this paper that the occurrence of a CFL 
quark matter core renders the compact star unstable. This confirms earlier 
findings by Baldo et al. [43J and Buballa et al. [44] for slightly different hybrid 
EsoS. 

An additional test to the mass-radius relationship is provided by a measure- 
ment of the gravitational redshift of line emissions from the CS surface. The 
disputed measurement of z = 0.35 for EXO 0748-676 |6J would be in accor- 
dance with both NMS and QCS interpretations. A measurement of z > 0.5 
could not be accommodated with the QCS model suggested here, while the 
NMS would not be invalidated by redshift measurements up to z = 0.6. 

Next we want to discuss the question whether measurements of the moment of 
inertia (Mol) / might serve as a tool to distinguish pure NMS from QCS mod- 
els. Due to the discovery of the relativistic double pulsar PSR J0737+3039 a 
measurement of the Mol became a possibility and has been recently discussed 
as another constraint on the EoS of compact stars, assuming that future mea- 
surements will exhibit an error of only about 10% [4"9~f50] . In our calculations 
we follow the definition of the Mol given in Ref. [TH] and show the results for 
the EsoS of the present paper in Fig. Due to the fact that the mass 1.338 
M Q of PSR J0737+3039 A is in the vicinity of the suggested critical mass 
region, the quark matter core is small and the expected Mol of the hybrid 
star will be practically indistinguishable from that of a pure hadronic one. 
The situation would improve if the Mol could be measured for more massive 
objects because the difference in the Mol of both alternative models for masses 
as high as 2 M Q could reach the 10% accuracy level. 

Finally we would like to discuss the question whether there are observables 
suited to distinguish between pure neutron stars and those with a quark matter 
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Fig. 2. Mass - radius relationship for CS sequences corresponding to a nuclear 
matter EoS (DBHF) and different hybrid star EsoS (DBHF+NJL), see text. Indi- 
cated are also the constraint on the mass from the pulsar J0751+1807 [2] and on 
the mass-radius relationship from the isolated neutron star RX J1856 [3] Present 
constraints on the mass-radius relation of CSs do not rule out hybrid stars. The dot- 
ted lines indicate the gravitational redshift, z = (1 — 2GM / R)~ l l 2 — 1, of photons 
emitted from the compact star surface. 

interior. As we have seen in the previous results, the hydrodynamic behavior of 
the hybrid star EoS has to be rather similar to that of pure hadronic matter 
in order to allow for sufficiently large masses. If so, the moment of inertia 
and other mechanical properties of the resulting stars will turn out to be 
indistinguishable to the level of a few percent. A different situation might 
occur for CS cooling where the transport properties and thus the excitation 
spectra of the dense matter play the essential role. As an example, pairing 
gaps of the order of an MeV or below will not affect the thermodynamics but 
are sufficiently large to influence on neutrino cooling processes. Let us discuss 
the example of the direct Urea (DU) process. 

If the DU process would occur in hadronic matter, it would give rise to a 
fast cooling and result in a strong sensitivity to slightest mass changes of the 
corresponding compact object. Therefore, it should not occur in CSs with 
masses below 1.5 M , as this would provide a cooling rate that is inconsistent 
with CS population syntheses [5T1I52] . If on the other hand the DU process 
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Fig. 3. Moment of inertia vs. the gravitational mass for the compact star EsoS 
discussed in the text. The highlighted mass regions correspond to the double pulsar 
J0737+3039 A+B and the pulsar J0751+1807. We also show anticipated data points 
with error bars corresponding to a measurement with a 1 a level of 10%. 

does not occur in a hadronic star, one would require that young, fast coolers 
such as Vela and 3C58 should have a rather large mass, again in contradiction 
with the present population syntheses [53] . 

A possible resolution to this hadronic cooling problem could be a phase transi- 
tion to quark matter with moderately enhanced cooling. This has been demon- 
strated for hybrid stars with a 2SC+X quark matter phase [54J which is in 
accordance with all presently known cooling constraints [52.53J. The physical 
nature of the hypothetical X-gap is, however, not yet clarified. A discussion 
of this issue can be found in Refs. [53JSB]. 

For the DBHF EoS the DU threshold is at n-^u = 0.375 fm -3 corresponding 
to a CS mass of M DU = 1.26 M Q , see Fig. HI The hybrid EoS presented in 
this work has a critical density for the quark matter phase transition which is 
below nrju provided a value rju > 1.024 is chosen. 

Thus for the parameter values r\y = 0.50 and rju > 1.024 the present EoS 
for hybrid star matter fulfills all modern observational constraints discussed 
above. 
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Fig. 4. Compact star masses vs. the central density for different values of the 
relative coupling strengths tjd and rjy, see Section [2J In a pure neutron (DBHF) 
star, DU processes occur for M > Mdu = 1-26 M (circle). The maximum mass 
of hybrid star configurations depends mainly on the relative strength of the vector 
meson coupling, rjy, while the transition density depends on rjy and rjr>- 



4 Summary 



We have investigated the compatibility of present constraints on the high 
density EoS with the concept of CSs possessing a QM core. The hadronic 
part of the EoS was taken from the DBHF approach, while the QM EoS is 
provided by a chiral NJL-type quark model with current- current interactions 
in the color singlet isoscalar scalar and vector meson channels as well as in 
the color antitriplet scalar diquark channels. The finite vector meson coupling 
enables us to describe large hybrid star masses by stiffening the QM EoS, 
whereas the chosen value for the diquark coupling ensures a sufficiently early 
phase transition to QM in order to avoid the activation of DU cooling within 
the hadronic shell of the hybrid star. Since the high density QM part of the EoS 
is softer than the corresponding pure hadronic EoS also the flow constraint is 
fulfilled in this scenario. We discussed the possibility to distinguish between 
pure NSs and hybrid stars by hypothetical data from successful measurements 
of the Mol for stars of a given mass. It turned out, that this would be possible 
for rather massive objects (M ~ 2.0 M ) provided the standard deviation of 
the measurements is less than 10%, as expected for PSR J0737+3039 A. 
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As our main result we conclude that no present phenomenological finding bears 
a strong argument against the presence of a QM core inside NSs. Moreover, 
we demonstrated that problems with cooling and flow which appear as weak 
points of a purely hadronic EoS at large densities can be resolved in a natural 
way when a transition to QM occurs at not too large densities. 

The most common argument against the presence of QM in CSs, resulting 
from the prejudice of the softness of QM EsoS, is no longer valid if we account 
for a vector meson interaction channel which stiffens the EoS. The earlier 
finding [43.25] that the occurrence of CFL quark matter renders hybrid stars 
unstable has been reconfirmed and therefore it should become customary for 
phenomenological scenarios of compact star evolution based on the assump- 
tion of a CFL quark matter phase (see, e.g., Refs. [39,40 41,42]) to check the 
necessary constraint of stability against gravitational collapse. 

The present study of dense quark matter in compact stars, albeit based on 
a microscopic dynamical approach used the strongly simplified NJL model 
interaction. It is desireable to repeat the present study for more realistic mi- 
croscopic approaches to the QCD EoS, e.g., on the basis of the QCD Dyson- 
Schwinger equation approach where interactions are nonlocal and momentum- 
dependent. One intermediate step towards this goal could consist in using 
Lattice QCD results on the momentum dependence of the quark selfenergies 
[5T] to adjust covariant formfactors of nonlocal, separable models [58.52] for 
the equation of state of quark matter in compact stars. 



Note added in proof 

After receiving the proofs for this Letter we have been informed that the mass 
of pulsar J0751*1807 has been corrected by Nice [HD] using a longer data base 
which allows to improve the value of the orbital decay and the detection of 
the Shapiro delay. Instead of previously 2.1 ± 0.2 M (la) the new value is 
1.26 M Q with estimated errors from 1.12 - 1.30 M (la) and 0.98 - 1.53 M 
(2<t). Therefore the mass of this object does no longer hold as a constraint 
on the EOS as described above. However, other new observations provide 
mounting evidence for the existence of a high maximum compact star mass. 
In addition to the LMXB 4U 1636-536 with 2.0 ± 0.1 M Q |4] mentioned in 
the introduction, very recent measurements on millisecond pulsars in globular 
clusters yield M = 1.96 + 0.09/ - 0.12 M and 2.73 ± 0.25 M (la) for PSR 
B1516+02B and PSR J1748-2021B, respectively [SI]. We also note that recent 
observations of the thermal polar cap emission of the quiescent LMXB X7 in 
the globular cluster 47 Tuc yield a rather large radius of 14.5 + 1.8/ — 1.6 km 
for a 1.4 solar mass compact star [62], which comes close to the lower bound 
on the mass-radius relation provided by RX J1856-3754 which is quoted in 
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our Letter. In summary we conclude that the large radii and masses found in 
several compact stars require a rather stiff equation of state. 



Acknowledgments 

We are grateful for critical reading and comments by M. Alford, M. Buballa, 
A. Drago, G. Grunfeld and F. Weber. We acknowledge the support by the 
Virtual Institute of the Helmholtz Association under grant VH-VI-041. The 
work of H.G. was supported in part by DFG grant No. 436 ARM 17/4/05 and 
that of F.S. by the Swedish Graduate School of Space Technology and the 
Royal Swedish Academy of Sciences. D.B. is supported by the Polish Ministry 
of Science and Higher Education. 



References 

[I] T. Klahn et al, Phys. Rev. C 74, 035802 (2006). 

[2] D.J. Nice, E.M. Splaver, I.H. Stairs, O. Lohmer, A. Jessner, M. Kramer, and 
J.M. Cordes, Astrophys. J. 634, 1242 (2005). 

[3] J. E. Triimper, V. Burwitz, F. Haberl and V. E. Zavlin, Nucl. Phys. Proc. Suppl. 
132, 560 (2004). 

[4] D. Barret, J. F. Olive and M. C. Miller, Mon. Not. Roy. Astron. Soc. 361, 855 
(2005). 

[5] F. Ozel, Nature 441, 1115 (2006). 

[6] J. Cottam, F. Paerels and M. Mendez, Nature 420, 51 (2002). 
[7] M. Mendez, private communication (2006). 

[8] M. Alford, D. Blaschke, A. Drago, T. Klahn, G. Pagliara and J. Schaffher- 
Bielich, Nature 445, E7 (2007). 

[9] M. Alford, M. Braby, M. W. Paris and S. Reddy, Astrophys. J. 629, 969 (2005). 

[10] A. Akmal, V.R. Pandharipande and D.G. Ravenhall, Phys. Rev. C 58, 1804 
(1998). 

[II] E.N.E. van Dalen, C. Fuchs, and A. Faessler, Nucl. Phys. A 744, 227 (2004); 
Phys. Rev. C 72, 065803 (2005). 

[12] E.N.E. van Dalen, C. Fuchs, and A. Faessler, Phys. Rev. Lett. 95, 022302 (2005). 

[13] T. Gross-Boelting, C. Fuchs, A. Faessler, Nucl. Phys. A 648, 105 (1999). 



13 



[14] P. Danielewicz, R. Lacey and W. G. Lynch, Science 298, 1592 (2002). 

[15] G. Stoicea et al. [FOPI Coll.], Phys. Rev. Lett. 92, 072303 (2004). 

[16] C. Fuchs, Prog. Part. Nucl. Phys. 56, 1 (2006); C. Fuchs et al., Phys. Rev. Lett. 
86, 1974 (2001); Ch. Hartnack, H. Oeschler, J. Aichelin, Phys. Rev. Lett. 96, 
012302 (2006). 

[17] D. Blaschke and D. Sedrakian (Eds.), "Superdense QCD Matter And Compact 
Stars", Dordrecht, Netherlands: Springer (2006). 

[18] N. K. Glendenning, "Compact stars: Nuclear physics, particle physics, and 
general relativity", New York, USA: Springer (2000). 

[19] F. Weber, "Pulsars as astrophysical laboratories for nuclear and particle 
physics," Bristol, UK: IOP (1999). 

[20] C. D. Roberts and A. G. Williams, Prog. Part. Nucl. Phys. 33, 477 (1994). 

[21] R. Alkofer and L. von Smekal, Phys. Rept. 353, 281 (2001). 

[22] D. Blaschke, C. D. Roberts and S. M. Schmidt, Phys. Lett. B 425, 232 (1998). 

[23] C. D. Roberts and S. M. Schmidt, Prog. Part. Nucl. Phys. 45, SI (2000). 

[24] D. Nickel, R. Alkofer and J. Wambach, Phys. Rev. D 74, 114015 (2006). 

[25] M. Buballa, Phys. Rept. 407 (2005) 205 

[26] M. Baldo, G. F. Burgio, P. Castorina, S. Plumari and D. Zappala, Phys. Rev. 
C 75, 035804 (2007). 

[27] S. M. Schmidt, D. Blaschke and Yu. L. Kalinovsky, Phys. Rev. C 50, 435 (1994). 

[28] D. Blaschke, G. Burau, M. K. Volkov and V. L. Yudichev, Phys. Atom. Nucl. 
62, 1919 (1999) [Yad. Fiz. 62, 2071 (1999)]. 

[29] S. Lawley, W. Bentz and A. W. Thomas, J. Phys. G 32, 667 (2006). 

[30] D. Blaschke, S. Fredriksson, H. Grigorian, A. M. Oztas and F. Sandin, Phys. 
Rev. D 72, 065020 (2005). 

[31] S. B. Riister, V. Werth, M. Buballa, I. A. Shovkovy and D. H. Rischke, Phys. 
Rev. D 72, 034004 (2005). 

[32] H. Abuki and T. Kunihiro, Nucl. Phys. A 768, 118 (2006). 

[33] H. J. Warringa, D. Boer and J. O. Andersen, Phys. Rev. D 72, 014015 (2005). 

[34] C. Gocke, D. Blaschke, A. Khalatyan and H. Grigorian, |arxiv:hep-ph /0104183. 

[35] F. Sandin and A. M. Oztas, Phys. Rev. C 73, 035203 (2006). 

[36] D. Blaschke, T. Klahn and D. N. Voskresensky, Astrophys. J. 533, 406 (2000). 

[37] D. Page, M. Prakash, J. M. Lattimer and A. Steiner, Phys. Rev. Lett. 85, 2048 
(2000). 



14 



[38] D. Blaschke, H. Grigorian and D. N. Voskresensky, Astron. Astrophys. 368, 561 
(2001). 

[39] M. Alford, P. Jotwani, C. Kouvaris, J. Kundu and K. Rajagopal, Phys. Rev. D 
71, 114011 (2005). 

[40] J. Berdermann, D. Blaschke, H. Grigorian and D. N. Voskresensky, Prog. Part. 
Nucl. Phys. 57, 334 (2006). 

[41] A. Drago, A. Lavagno and G. Pagliara, Nucl. Phys. A 774, 823 (2006). 

[42] D. Page and A. Cumming, Astrophys. J. 635, L157 (2005). 

[43] M. Baldo, M. Buballa, F. Burgio, F. Neumann, M. Oertel and H. J. Schulze, 
Phys. Lett. B 562, 153 (2003). 

[44] M. Buballa, F. Neumann, M. Oertel and I. Shovkovy, Phys. Lett. B 595, 36 

(2004) . 

[45] H. Grigorian, Phys. Part. Nucl. Lett. 4, 382 (2007). 

[46] G. R6pke, D. Blaschke and H. Schulz, Phys. Rev. D 34, 3499 (1986). 

[47] D. Blaschke, T. Tovmasian and B. Kampfer, Sov. J. Nucl. Phys. 52, 675 (1990); 
[Yad. Fiz. 52, 1059 (1990)]. 

[48] D. N. Voskresensky, M. Yasuhira and T. Tatsumi, Nucl. Phys. A 723, 291 
(2003). 

[49] M. Bejger, T. Bulik and P. Haensel, Mon. Not. Roy. Astron. Soc. 364, 635 

(2005) . 

[50] J. M. Lattimer and B. F. Schutz, Astrophys. J. 629, 979 (2005). 

[51] S. Popov, H. Grigorian, R. Turolla and D. Blaschke, A&A 448, 327 (2006). 

[52] S. Popov, H. Grigorian and D. Blaschke, Phys. Rev. C 74, 025803 (2006). 

[53] D. Blaschke and H. Grigorian, Prog. Part. Nucl. Phys. 59, 139 (2007). 

[54] H. Grigorian, D. Blaschke and D. Voskresensky, Phys. Rev. C 71, 045801 (2005). 

[55] D. Blaschke, D. N. Voskresensky and H. Grigorian, Nucl. Phys. A 774, 815 

(2006) . 

[56] D. Blaschke, PoS JHW2005, 003 (2006). 

[57] M. B. Parappilly, P. O. Bowman, U. M. Heller, D. B. Leinweber, A. G. Williams 
and J. B. Zhang, Phys. Rev. D 73, 054504 (2006). 

[58] D. Blaschke, H. Grigorian, A. Khalatyan and D. N. Voskresensky, Nucl. Phys. 
Proc. Suppl. 141, 137 (2005). 

[59] D. Gomez Dumm, D. B. Blaschke, A. G. Grunfeld and N. N. Scoccola, Phys. 
Rev. D 73, 114019 (2006). 



15 



[60] Conference contribution by D. Nice, 40 Years of Pulsars, Montreal, 2007. 

[61] Conference contribution by P.C.C. Freire, et al., 40 Years of Pulsars, Montreal, 
2007. 

[62] G. B. Rybicki, C. O. Heinke, R. Narayan and J. E. Grindlay, Astrophys. J. 644, 
1090 (2006). 



16 



